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We derive inclusive and diffractive forward-gluon production in the scattering of a qq dipole off an 
arbitrary target in the high-energy eikonal approximation, suitable to study the saturation regime. 
We show how the inclusive cross-section is related to the total cross-section for the scattering of 
a colorless pair of gluons on the target: the gluon-production cross-section can be expressed as a 
convolution between this gg dipole total cross-section and a dipole distribution. We then consider as 
an application the forward-jet production from an incident hadron and describe forward-jet produc- 
tion at HERA and Mueller-Navelet jets at Tevatron or LHC. We show how these measurements are 
related to the qq — gg or gg — gg dipole-dipole cross-sections and why they are therefore well-suited 
for studying high-energy scattering in QCD. 



I. INTRODUCTION 

o ■ 

Understanding high-energy scattering in QCD has been the purpose of a lot of work over the past ten years. At 
£S| [ high energies, as one approaches the unitarity limit, high-densityphases of partons are created and non-linear effects 
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become important. In this regime called saturation P, S S ; the basic properties of perturbative QCD, such 
as factorization and linear evolution, breakdown. To study scattering near the unitarity limit, the qq dipole model 
0,0] has been developed. This formalism constructs the light-cone wavefunction of a dipole (a quark-antiquark pair 
in the color singlet state) in the leading logarithmic approximation. The size of the dipole provides a scale which 
is supposed to be much smaller than 1/Aqcd to justify the use of perturbation theory. As the energy increases, 
the dipole evolves and the wavefunction of this evolved dipole is described as a system of elementary dipoles. This 
formalism is well-suited to study high-energy scattering because when this system of dipoles scatters on a target, 
density effects and non-linearities that lead to saturation and unitarization of the scattering amplitude can be taken 
^T) | hrto account. 

On a phenomenological point of view, the key point is to relate the experimental probes to the scattering of a dipole. 
In the case of lepton-hadron collisions this is straightforward: a lepton undergoes hadronic interactions via a virtual 
photon and one can view the interation as the splitting of the virtual photon into a qq pair which then interacts. 
This dipole picture has had great succes for the phenomenology of hard processes initiated by virtual photons. In 
hadron-hadron collisions however, there are no such virtual-photon probes. The probes are the final-state particules 
that we measure, e.g. a final-state gluon that we detect as a jet. Establishing a link between the scattering of a 
dipole and observables like jet cross-sections is of great theoretical and phenomenological interest in the prospect of 
the LHC, if one considers the impact that the dipole picture has had on HERA phenomenology for example and how 
O ' the HERA measurements have helped understanding high-energy scattering. 

First steps in this direction were taken in [T^j where gluon radiation from a quark was calculated. In , Kovchegov 
' and Tuchin derived the cross-section for inclusive gluon production in deep inelastic scaterring off a large nucleus. In 
both cases, the cross-sections are expressed in terms of the scattering of a gg dipole (a gluon-gluon pair in the color 
singlet state) on the target. It is crucial in jTlj that the target is a nucleus since they work in an approximation 
in which all the scatterings on the nucleons happen via one or two-gluon exchanges |l2j. On a phenomenological 



side, descriptions of Mueller-Navelet jets in 13j and forward jets in [lj| are performed in terms of qq dipole-qq 
dipole scattering, leading to interesting results and predictions. The key point in these analysis is the description of 
the jet vertex in terms of qq dipoles: it is done using the equivalence between the kr— factorization and the dipole 
factorization. The factorization formula used is therefore not justified when the dipole-dipole cross-section is not 
kx— factorizablc which could be the case at high energies. 

In this paper, we consider the collision between a qq dipole and an arbitrary target and calculate the cross-section 
for the production of a gluon in the forward direction of the dipole. We work in the eikonal approximation valid 
for high-energy scattering in which the interaction with the target is described by Wilson lines. We derive formulae 
for the inclusive and diffractive cross-sections. We express the inclusive cross-section in terms of the scattering of a 
colorless gluon-gluon pair on the target. We generalize the result of [ 1 1| . obtaining it in a very different framework and 
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for the scattering off any target. We write the inclusive cross-section in the factorized form given below in equation 
(|40|l : a convolution between an effective dipole distribution and the total cross-section for the scattering of a gg dipole 
on the target; such a factorization does not happen for the diffractive cross-section. We obtain a very general formula 
that includes all numbers of gluon exchanges and non-linear quantum evolution. 

We then take another step towards the description of observables in hadronic collisions as we extend our results to 
the case of an incident hadron instead of an incident qq pair. This is done by including linear evolution before the 
emission of the measured gluon and by using collinear factorization. We thus describe forward-jet cross-sections at 
hadronic colliders in terms of the scattering of a gg dipole on the target, see later formula (|55|l which exhibits a dipole 
factorization for the forward-jet probe and generalizes the formulation of 0,0]- 

Finally, as an application to our formulae, we consider the target to be a virtual photon, or equivalently a qq dipole, 
and obtain the forward-jet cross-section at HERA |Tfij . Its high-energy behavior is driven by the quantum evolution 
of the gg dipole-qq dipole scattering, making the forward-jet measurement of particular interest to investigate high- 
energy scattering and unitarization in QCD. The same conclusion holds for the Mueller-Navelet jet cross-section [T(| 
at Tevatron or LHC that we also derive. 

The plan of the paper is as follows. In Section 2, single inclusive forward- gluon production in the scattering of 
a qq dipole off an arbitrary target is calculated. The diffractive cross-section is derived in Section 3. The relation 
between the inclusive cross-section and the gg dipole-target scattering is derived in Section 4. In section 5, extending 
the previous result to the case of an incident hadron and specifying different targets, we derive the Mueller-Navelet 
jet and forward-jet cross-sections. Section 6 concludes. 



II. INCLUSIVE GLUON PRODUCTION 




FIG. 1: Inclusive gluon production off a qq dipole. xo and x\ : transverse coordinates of the quark and antiquark; k and y : 
transverse momentum and rapidity of the measured gluon. The black points represent emission of the gluon by either the quark 
or the antiquark. The vertical wavy lines represents the interaction with the target. The gluon is emitted before the interaction 
or after the interaction in which case it comes with a minus sign as explained in the text. 



In this section we derive single inclusive gluon production in the high-energy scattering of a qq dipole off an arbitrary 
target. We shall use light-cone coordinates with the incoming dipole being a right mover and work in the light-cone 
gauge A + = 0. In this case, when the dipole passes through the target and interacts with its gauge fields, the 
dominant couplings are eikonal: the partonic components of the dipole (a quark, an antiquark and soft gluons) have 
frozen transverse coordinates and the gluon fields of the target do not vary during the interaction. This is justified 
since the incident dipole propagate at nearly the speed of light and its time of propagation through the target is 
shorter than the natural time scale on which the target fields vary. The effect of the interaction with the target is 
that the components of the dipole wavefunction pick up eikonal phases. 

Let us be more specific. In Fig.l is represented the production of a gluon of transverse momentum k and rapidity 
y off a dipole with transverse coordinates xq for the quark and X\ for the antiquark. The size of the dipole \x$ — x\\ 
is supposed to be small in order to justify the use of perturbation theory (|a;o~ x\\ <C I/Aqcd)- We work in a 
frame in which the dipole rapidity is not too large so that the radiation of extra softer gluons is to be described by 



quantum evolution of the target. The necessity of the last term can be seen by considering the case where there are 
no interactions: the emission-after-interaction term is necessary to get zero gluon production without interaction. 
The incident hadronic state is a colorless qq pair and has the following decomposition on the Fock states: 



\d) = \d) + \dg) • 



The bare dipole \d)o is characterised by the wavefunction 



(i) 



(2) 



where a and a denote the colors of the quark and antiquark respectively and where the transverse positions of the 
partons have been specified. The qqg part of the dressed dipole \dg)o is characterised by the wavefunction 



\dg)o 
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where the gluon is characterized by its color a, its polarization A, its rapidity y and its transverse coordinate z. e\ is 
the transverse component of the gluon polarization vector and T a is a generator of the fundamental representation of 
SU(N C ). The term in brackets in @ is the well-known wavefunction for the emission of a gluon off a qq dipole 9]; 
the two contributions correspond to emission by the quark and antiquark respectively. The only assumption made 
to write down Q is that the gluon is soft, that is its longitudinal fraction of momentum with respect to the incident 
dipole is small. As already mentionned, we work in a frame in which only bare or one-gluon components need to be 
considered in the wavefunction \d); softer gluons will be included through the quantum evolution of the target. One 
then truncates the perturbative expansion Q at first order in the strong coupling constant g s . 

Let us denote the initial state of the target \t). The outgoing state is obtained from the incoming state \ d)<S ) \t) by 
action of the S— matrix. In the eikonal approximation, S acts on quarks and gluons as (see for example 17, 18, 19]): 



S\(a,x))®\t) = [W F {x)] aa , | (a',*))® |*) , 



S\(a,\,z,y))®\t) =Y,W% b (z)\(b,\,z,y))®\t) , (4) 



where the phase shifts due to the interaction are described by Wf and Wa, the eikonal Wilson lines in the fundamental 
and adjoint representations respectively, corresponding to propagating quarks and gluons. They are given by 



W F , A {x) = V exp{ig s J dz+T^ A A a _(x, z+)} 



(5) 



with A- the gauge field of the target and Tp A the generators of SU(N C ) in the fundamental (F) or adjoint (A) 
representations. V denotes an ordering in z + . 

Therefore the state \^ ut) =S\d)®\t) emerging from the eikonal interaction reads \^ ou t) = l^i) + 1 1 ^) with the 
componants given by 
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|^2) represents the first contribution pictured in Fig.l while the second contribution is hidden in To see that, 

we must express the Fock states which appear in |^ ou t) in terms of physical states. At first order in g s one has: 



\(a,x ),(a,xi)) = \(a,x ),(a,Xi))phy S -'^2 J ^zdy 

fiXa 
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\(a,Xo), (a, xi), (b, A, z,y)) = \(a,x ),(a,x 1 ),(b 7 X,z,y)) phys 



(9) 



One immediately sees that the emission-after-interaction term arises with a minus sign. Dropping the qq part which 
does not conribute to gluon production, using ©, 0, ©, and ©, the out going state can be rewritten 
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The different contributions contained in this wavefunction have a straightforward physical meaning: the term contain- 
ing e\.(z— xq) comes from the emission of the gluon by the quark while the one containing e\.{z— x\) comes from to 
the emission by the antiquark. Moreover for both these terms, the contribution with three Wilson lines corresponds 
to the interation happening after the gluon emission while the contribution with two Wilson lines corresponds to the 
interation happening before. 

From the outgoing state 1)10(1 , the gluon-production cross-section reads 
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where a\ x (k,y) and a Ci \(k,y) are respectively the creation and annihilation operators of a gluon with color c, polar- 
ization A, rapidity y and transverse momentum k. xoi = x$— x\ is the size of the incoming dipole and b— {xq+X\)/2 is 
the impact parameter. Let us rewrite the cross-section 1)111 using operators a c ,x(z,y) in transverse coordinate space 
that act on \^ ut) '■ 
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d 2 b e ik -^{t\P(x ,x 1 ,z 1 ,z 2 )\t) 



Z\ and z 2 represent now the transverse coordinates of the measured gluon in the amplitude and the complex conjugate 
amplitude respectively. P(xq, Xi,Z\,Z2) can be easily calculated from 1)10)1 and its definition 1)12)1 using 



a c ,\( z > 2/)K a ' ^o), (a, (c', A', z', y')) = S cc >8\x/5 ( - 2 " 1 (z-z')5(y-y')\(a, x ), (a, xx)) 
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It is given by 

P{x ,xi,z 1 ,z 2 ) = Tr 
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This expression contains terms with four, five, or six Wilson lines, however it can be reduced quite easily using the 
following formulae: 



Tr(wl{x)T a W F {x)T b ^Wf{z) = ^Tr (w A (x)W A (z) 



(15) 



Tr 



{w F {x l )T a wl{x x )W F {x )T a wl{x )) = ^Tt{w\( Xi )W a (xo)) . (16) 



These are obtained using the Fierz identity that relates the Wilson lines Wf and Wa '■ 

[W F (x)] tJ [wUx)]ki = ^ 6 a S jh + 2W?(x)T? l T* j . (17) 

In formula (|14|l . terms which a priori contain six Wilson lines reduce to traces of two adjoint Wilson lines as 
W f (x)Wf(x) — 1. For the same reason, terms which display five Wilson lines contain only three and can be re- 
duced using l|15(l . Finally, terms with four Wilson lines arc either trivial and equal to CfN c or can be reduced using 
(|16fl . At the end, the cross-section (I12fl obtained from P(xq ,Xi,Zi,Z2) contains only traces of two adjoint Wilson lines 
and is given by 

Tr^^^A^i)))^^^^)^^^))^^^^)^^)))^ (18) 



cPkdy 



where we have denoted (t\ . \t) = ( . )f. 

Interestingly enough, the only quantity involved in l|18|) is the forward scattering amplitude of a colorless pair of 
gluons on the target: 

T(x,x\Y-y) = 1- j^(Tr{w{(x')W A (x))) t . (19) 

This scattering amplitude is present in (|18fl for all the dipoles involved in the process. It contains the scatterings 
with all numbers of gluon exchanges and, via its quantum evolution, it also contains the emissions of gluons softer 
than the measured one (k,y). Indeed, if Y is the total rapidity, the amplitude (fH^l depends on Y — y. The cross- 
section <|18[) is made of four contributions, each of them proportional to one of the (a;, — zi).(xj — z^) (i,j — 0, 1), for 
which the gluon is emitted from the (anti) quark at transverse coordinate Xi in the amplitude and Xj in the complex 
conjugate amplitude. Each contribution itself contains four amplitudes T whose physical origins are the following: 
T(xi, Xj, Y—y) corresponds to the gluon being emitted after the interaction both in the amplitude and in the complex 
conjugate amplitude in which case it does not interact; T{xi, Z2, Y—y) corresponds to emissions after the interaction in 
the amplitude and before the interaction in the complex conjugate amplitude and vice- versa for T(z\, Xj, Y—y); finally 
T(z±, Z2, Y — y) corresponds to emissions before the interaction both in the amplitude and in the complex conjugate 
amplitude. 

The scattering amplitude T{x, x' , Y—y) is related to the total cross-section ai gg \ t [x—x' , Y—y) for the scattering of 
a gluon dipole (gg) of size x — x' on the target, with rapidity Y — y: 

<j {gg)t (x-x',Y-y) = 2 J d 2 [^pjT{x,x',Y-y) . (20) 

To compute this cross-section, one has to evaluate the target averaging ( . ) t contained in T (see formula (|19ll ) which 
amounts to calculating averages of Wilson lines in the target wavefunction. A lot of studies are devoted to this 
problem as we briefly discuss later in Section V. Here we only establish the link between the observable l|18fl and the 
dipole amplitude ltT9"|) . 

Formula l|18l) generalizes the result of ITll l20j where the same cross-section has been derived for a target nucleus 
in the quasi-classical approximation of |l2| in which the scatterings on each nucleon happen via one or two gluon 
exchanges. It is quite remarkable that only one dynamical quantity appears in the result and this is what will allow 
us to write the inclusive cross-section l|18t in a factorized form as is shown in Section 4. But first, we shall calculate 
the diffractive gluon-production cross-section to exhibit the difference with the inclusive one and see that in this case, 
more than one dynamical quantity is involved. 



III. DIFFRACTIVE GLUON PRODUCTION 



The diffractive cross-section for gluon production in a qq dipole-target collision is pictured in Fig. 2 where we define 
the diffractive process as one in which the outgoing wavefunction is in a color singlet state and in which the target 
does not break up. The cross-section is calculated in the same way than the inclusive one in the previous section. The 
difference is that we now have to project the outgoing state \^ ou t) on the subspace of color-singlet states. The final 
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FIG. 2: Difltactive gluon production off a qq dipole. The object exchanged between the qq pair and the target is a color singlet 
and the target does not break up. Notations: see Fig.l. 



state containing diffractive gluon production is \^ s ing) = PsingY^ out) where the projector on the color-singlet states 
P sing is given by: 

Psing = jt- I d 2 x d 2 x'( ^S^K/x, x), (u,x')) j ( § a , & ((a, x), (a,x')\ J + 1 ^ I ' d 2 x d 2 x' d 2 zdy 

J2T^,x),{u,x'),{b,X,z,y)) ( T^((a,x),(a,x%(a,\,z,y)\\ . (21) 

u,v,b I \a,a,a / 



One obtains with 0, 0, and (JSTJ: 



P S ing\Vi) = ^Tr^Wlix^Wpixo))^® \t) , 



(22) 
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Tr(wl(xi)T b W F (x )T c )w b A c (z) 



T£ a \(a,xo),(a,xi)Aa,\z,y))(g>\t) . (23) 



Psing 1^2) is the contribution of the emission-before-interaction term (first graph in Fig. 2). Writing \d)o = \d)— \dg)o in 
Psing makes the emission-after-interaction term arise (second graph in Fig. 2). One is then able to express losing) 
in terms of physical states. Keeping only the qqg part, one has: 



w = 4f £ / ^ v 



W, 
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TT 



e\.{xp-z) _ e x .{xi-z) 
\x -z\ 2 



\Xl — Z[ 



C 



^Tr(wU Xl )T b W F (x )T c )w h A c (z) 



-^ c ^( W k^)WF(x ))^T^ a \(a,x ),(a,x 1 ),{a,X,z,y)) phvs ^ \t) . (24) 

losing) is not exactly the final state one needs to consider since our definition of diffractive implies that the target 
does not break up and we have imposed that condition yet. The final state lf2*4^l would be the one to use to calculate 
diffractive gluon production without any requirement on the target. To insure the target not to break up, one 
projects the outgoing state losing) on the subspace spanned by \t). The effect is that the Wilson lines in become 
target-averaged: 



1 
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with 



$(z) = (Tr(wl(: Xl )T a W F {x a )T b )w% b (z)) t - C F (Ti(w f { Xi )W f {x )) ^ . (26) 

The diffractive gluon-production cross-section is now obtained from formula with Y^diff) instead of |^ 0U {). 
The final result is 



d 2 /fcdy v u ' ir 2 C F N 2 7 7 2tt 2vr 
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(27) 

Note that if one had computed the cross-section with losing), there would have been a global target average at the 
cross-section level and not averages at the amplitude level as is the case in (I27|l . The first term in (|26f) is proportional 
to the elastic S— matrix for the scattering of a colorless qqg triplet on the target while the second term is proportional 
to the elastic S— matrix for the scattering of a qq dipole. There is then two dynamical quantities playing a role in 
this diffractive cross-section. This is the main difference compared to the inclusive case in which there was only one. 

If one considers that the target is a nucleus, and that each scattering on the nucleons happens via a two-gluon 
exchange, then the target averages in l|26ll are computable (see e.g. 18]) and one recovers the result of 20]. Formulae 
()27J) and Ij26(l are a generalization to an arbitrary target including all numbers of gluon exchanges. 

Let us mention also that, making use of the following identity: 

2Ti(wl(x 1 )T a W F {x Q )T b ^Wt(z) = Tr(wUxi)W F (z)^Tr(w F (z)W F (x Q )) - (a>i)W>.(a;o)) , (28) 

one is able to relate the first term of H26(l to the scattering of two dipoles and obtain: 

$(z) = ~ (Tr(wl(x 1 )W F (z)^Tr(w F (z)W F (x ))) t - ^ (tv(wI( Xi )W f (x )^ ^ . (29) 

The first term in $ is now proportional to the elastic 5— matrix for the scattering of two qq dipoles on the target. 
What is usualy considered is the scattering off a large nucleus in which case one can justify neglecting the fluctuations 
and writing 

Tr(w F ( Xl )W F (z^TY(wl(z)W F (x )^ t = (Tt (w f { Xi )W f (z)) } (Tr (wl(z)W F (x j) ) . (30) 

Making this simplification is very useful to do phenomenology pi) , as it enables to deal with only one quantity but it 
is not correct for any target. Note also that if one writes the product in the ri ght- hand side of (|30() using T amplitudes 
as in H19(l , one recovers the two-gluon exchange approximation calculated in |22L l23| by neglecting the T 2 term. 



IV. DIPOLE- FACTORIZED FORM OF THE INCLUSIVE CROSS-SECTION 

We now come back to the case of the inclusive cross-section Ijl8|l and show in this section that it can be significantly 
simplified and expressed in a form that exhibits a dipole factorization. 

On the right-hand side of 118fl . the term in brakets contains four amplitudes and for each of them two of the 
integrations can easily be carried out. Let us begin with the amplitude which does not depend on 6, the integration 
of the emission kernel over that impact parameter gives: 

d ° I 121 12 = 27Tlo S] 1 1 (31) 

\x i -zi\ z \x j -z 2 \' ! \ Xl j-Zl + Z 2 \ 

where and p is an infrared cutoff which will eventually disappear. Then one writes: 

e >H. nh ^_JL_^y. v) _ I J £ e—log^ „ M ,(z,Y- v ) . (32, 



Consider now the amplitude which only depends on z\, we then perform the z 2 integration 

d 2 Z2 C tk.(z 2 - Zl ) {x i -z 1 ).{x j -z 2 ) = k.(xj- Zl ) j k . {xj - Zl) ^ 
27r \xi-zi\ 2 \xj-z 2 \ 2 ik 2 \xi-zi\ 2 



and obtain 

One can easily see that the amplitude which only depend on z 2 gives the same result. For the last amplitude which 
does not depend either on z\ or Z2, we integrate the emission kernel: 

d 2 z x d 2 z 2 jk.( Z2 - Zl ) {x i -z 1 ).(x j -z 2 ) = 1 ik , Xi . . , 

2ir 2tt \x. l -z 1 \ 2 \x j -z 2 \ 2 k 2 ( ' 

and the b integration of the amplitude yields to the cross-section. Putting the pieces together, the gluon production 
cross-section reads 

Jk M ■-** ± <-*r>(**f*- - 1 *»to-*)^,r-v) . (3e, 

Performing the summation over i and j, the cutoff p disappears and the two terms displaying 8^ (z) also vanish 
because <J( flfl )t(0, y) = 0. Using then the following results for two dimensionnal vectors where 6 is the Heavyside step 
function: 

. 2 \e x log- ^— =47re(M-|z|) log M , (37) 

o \x — z\\x + z\ \x\ 



2 " 2i f n k.z k.(z-x) k.(z+x)\ 8in „ , , .. 

^ " itiF - - |^J = m cob WN-W) , (38) 

one carries out the integration over the azimutal angle of xoi- Denoting the size of the incident qq dipole \xqi \ =tq, 
the result reads 

^ M -^J*'^"{(^~>-*fi)^-W) + p<*-tf^Y-v). (39) 

Most of the time the cross-section cri gg \ t {z,y) will not depend on 9 as one averages over the azimutal angle of the 
target. This enables the final 9 integration: 

do , . 1ol s Cf f n / j\ , , , , 

( r °) = uT~ I dz nro,z,k)o igg)t (z,Y-y) . (40) 
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We have singled out the following normalized distribution 



"<r„. :./•■) Hlr„ :\[- z J 1 {kz)-^J {kz)\og{r /z))+J Q {kz)6(ii- ^,. Ml) 



dz 2 <p(ro, z, k) = 1 



(42) 



Note that now, z and k respectively stand for \z\ and \k\. The gluon production cross-section as written in (|4C)|> exhibits 
a dipole factorization: the gg dipole-target cross-section (J( gg )t is convoluted with the effective dipole distribution cj) 
which contains the emission of the forward jet of momentum k off the initial qq dipole of size rg. It is remarkable that 
we obtain factorization and such a simple formula. 
One can easily check that <f> (j41|) verifies: 



2z4>(r , z, k) = 



d 



6(r - z)z-^- (jo(kz) log — 
oz\ z 



Thanks to this identity, one can express the gluon production cross-section (|40|) in another convenient way: 



da 
d 2 kdy 



M 



2a s C F p r d ( d 

— J o d zJ {kz)lo g --[z-o {gg)t {z,Y-y) 



(43) 



(44) 



Formulae l|40(l and 144(1 express the single inclusive forward-gluon production cross-section in the scattering of a qq 
dipole off an arbitrary target and are the main results of the paper. The intermediate cross-section crr gg \ t for the 
scattering of a gg dipole on the target includes all multiple scatterings and linear or non-linear quantum evolution. 
As shown by (|19l) and (|2Up. it can be computed from a trace of Wilson lines correctly averaged over the target gluon 
fields. 

So far, we have only considered the case where the measured gluon is the closest in rapidity to the dijet coming 
from the dissociation of the incident qq pair. One can extend the result to the case of a gluon measured at a less 
forward rapidity by including the emissions of harder gluons, emitted before the gluon detected in the final state. 
To stay consistent with our soft-gluon approximation, this can be done by including leading logarithmic evolution 
before the emission of the measured gluon. This is quite straightforward if one sticks to linear BFKL evolution which 
is justified as long as the rapidity of the measured gluon is not to big and non-linearities are negligible. Since the 
incident particule is a qq dipole, this BFKL evolution is the usual dipole evolution from the initial dipole of size r to 
the dipole of size r from which the measured gluon is now emitted. In formulae, that amounts to do the substitution: 

Jhr y [ro) ^ I ^ n{ro ' r ' y) Jk {r) (45) 

where n(ro, r, J/Ljs the number density of dipoles of size r inside the initial dipole of size r$ at rapidity y. This quantity 
was defined in and was shown to satisfy the BFKL equation [24|. The solution is: 

n(r , r ' y ^ = 2 J ^ ("") exp (~~yx(7)j (46) 

where the complex integral runs along the imaginary axis from 1/2 — zoo to 1/2 + ioo and where the BFKL kernel is 
^(7) = 2^)(1) — ^(7) — ^(1— 7). The substitution (gHJ gives: 

ffijfijW = J dz Mkz)n(r ,z,y)-^- a (gg)t (z,Y-y) j (47) 

where 

The most natural application to our formulae is a measurement in deep inelastic scattering in which one detects 
the most closest jet to the dijet coming from the dissociation of the photon in case of formula 144f> or a less forward 
jet for formula (|47|l . Such a measurements would certainly give information on the dipole- proton scattering. However 
the rest of our study will not concentrate on that, it will turn to the case of forward jets emitted by hadrons. Such 
jets represent the most natural probes in hadronic collisions and are a very interesting application to our formulae. 

V. AN APPLICATION: FORWARD-JET PRODUCTION 

In this part we consider the same cross-section as in the previous section, but with the initial particule being a 
hadron. Indeed in experiments, collisions that procuce jets are proton-lepton or proton-antiproton collisions. Forward 
jets are jets measured in the forward directions of the hadrons. We cannot compute such processes directly because 
of the non-perturbative nature of the hadrons. Therefore we start by explaining how to extend the previous result 
from an incident perturbative qq dipole to an incoming hadron. In a second part, we also specify the target and 
consider a virtual photon, the goal being to describe forward jets in deep inelastic scaterring. We finally generalize to 
Mucller-Navclct jets at hadron colliders. 

A. Gluon production off a hadron 

In this section we explain how to make the link between the qq dipole and the hadron using collincar factorization. 
Indeed, providing that the transverse momentum of the measured jet is larger than a perturbative scale, forward-jet 
production is a hard process and obeys collinear factorization. Moreover, since we are dealing with forward jets, the 
final state gluon has a fraction of longitudinal momentum (with respect to the incident hadron) xj that is never too 
low, thus there are no small— a; effects in the wavefunction of the hadron that could break collinear factorization. 



Here is how we proceed. Consider first the cross-section (1471 derived in the previous section and take the collinear 
limit for the incident dipole tq^l 3> 1 where fi is an arbitrary hard scale called the factorization scale. This limit 
requires the emission of the gluon to occur at a scale k~ \i much larger than the scale of the initial dipole l/Vo- Such 
a strong ordering between the two scales is indeed what leads to the collinear factorization. Moreover it is justified in 
this case because we really want to describe the emission of the jet off a hadron which naturally carries a soft scale. 
In this appropriate limit, h(ro,z, y) reduces to: 

-/ n f d l ( r oM) 27 (2a s C F \ 2 f dj 2 I (2a s C F \ 

since for r$n^>l the integral is determined by the 7~0 behavior of the integrand. We have used x{l) — 1/7 for 7 — 0- 
Then one recognizes that the right-hand side of l|49f) is an integral representation of a modified Bessel function [25j : 



dl u iA 2 7 1 _ ( 2a s C F \ _ 1 / 7rlog(rg^ 2 ) T ( n \ 2a s C F 1 , 



2 i »<-»"» 'V exp l^r 9 J - W^r /,|2,/ - log< '»' ,2)! '^ (50) 



And (|47|l reduces finally to: 




2a s C F log(r^ 2 ) r ( /2a s C F ,,.>., . 



d f d 

dz J (fcz) — f z— (T (s9)t (^, A77) ) (51) 



where we have changed the variable y = \og\/xj and defined Arj — Y^ log l/xj, the rapidity interval between the jet 
and the target. The factor in brackets in (|51l) is the gluon distribution function xjgd(xj, fi 2 ) inside the incident qq 
dipole of size tq at the Double Leading Logarithmic (DLL) approximation. This is not surprising: we have assumed 
that the emitted gluon was soft which is responsible for the leading \og(l/x,j) and we have considered the collinear 
limit which is responsible for the leading log(rQ/i 2 ). Within this DLL approximation, the distribution gd inside the 
qq pair of size vq is the same than the gluon distribution inside a gluon of virtuality 1/r 2 with 2 C F being the color 
factor instead of Ca ■ Now the link with the gluon production cross-section off an incident hadron is obvious: the hard 
part of the cross-section is the same and the hadron appears in the cross-section via its parton distribution function 
One has then 

da h+t ^ J+x 1 f d ( d \ 

gh{xj^ 2 ) \ dz J (kz) — hs— <T(gg) t (z, At})) . (52) 



d 2 kdxj 2k 2 'J — "u^~> dz \ d 

Note that this formula is only valid at the DLL approximation since we have always taken the emitted gluon to be 
soft, that is a;jCl. 

A simple way to see that ()52(l is the correct formula is to consider the case of linear leading-logarithmic BFKL 
evolution between the jet and the target. In this case &( gg )t obeys /cy— factorization and can be expressed in terms of 
the unintegrated distribution in the target f F ■ 

af™\z,Ar ) ) = ^± j ^(l-Mqz))f T ( q ,Q Tl A v ) (53) 

where Qt is a scale characterising the target. Inserting Ij53|l into l|52(l gives 

lh+t-^J+X A AT 

9|g^ = ^ 9h ( X j,, 2 )Mk,Q T ,Arj) (54) 

and we recover the known BKFL result 26] for forward-jet production. 

Formula (|52|) can be written in a very simple factorized form if one considers the cross-section integrated over the 
transverse momentum of the jet above a given cut Et ■ 

dx~j ^ J 2 rffc2 d?kdxj =5/l ^ J ' /i2 ) J d%z $j( z i E T)°{ sg )t{zi A??) , (55) 

where the distribution <f>j is related to a bessel function: 



<f>j{z,E T ) = -Lj^Erz). 
2nz 



(56) 



We have assumed that the gg dipole-target cross-section <?( gg )t{z) goes to a constant as z goes to infinity fast enough 
to cancel the edge terms in the integrations by parts. Typically, it works with vtgg)t(z) ~ 1 — exp(— z 2 ) or with 
a (gg)t( z ) ~ 1 — exp (— log 2 (z)). Formula i(55() exhibits the dipole factorization for a forward jet emitted off a hadron. 
We recover the result obtained from indirect calculations in using the equivalence between the kr— factorization 
and the dipole factorization. The feature that this approach couldn't reveal in is that it is a gg dipole that has to 
be considered and not a qq dipole. We have also shown here that the dipole factorization is still valid when the 
cross-section G( gg )t is not kx— factorizable, justifying the use of 1(55(1 to study saturation effects in forward jets |l3lll4| . 



B. Forward jets in deep inelastic scattering 

The study of forward jets in deep inelastic scattering represents the most straightforward application of our formula. 
In such processes, a proton scatters on an electron and a jet is detected in the forward direction of the proton. Therefore 
the interaction takes place between the proton and a transversally (T) or longitudinally (X) polarized virtual photon, 
see Fig. 3a. An hadronic interaction with a virtual photon is easily expressed in terms of qq dipoles: the photon 
fluctuates into the quark-antiquark pair which then interacts. The wave functions <pr and <pL describing the splitting 
are well-known. The cross-section is then given by formula (|55|l with the gg dipo\e-qq dipole cross-section <J( gg )d 
carrying the rapidity dependance and with an extra convolution with 4>t,l ■ 

^ = 9p(xj,fi ) J d rd z 4> T x(r,Q )cj>j(z, E T )a {gg)d (r, z, Arj) . (57) 

The virtuality of the photon has been denoted Q 2 and Et is the cut on the jet transverse momentum. In the regime 
where Q~£'t, DGLAP evolution is suppressed and fixed-scale quantum evolution drives the high-energy behavior of 
the cross-section. This evolution is contained in &( gg )d which is the only unknown quantity in formula 157|) . It can 
be obtained by computing a trace of Wilson lines averaged over the fields of the dipole, see formulae ((19(1 and l|2()|l: 
performing the target averages and calculating ct( gg )d is however not easy. The color glass condensate theory provides 
a way of doing it: 

(Tr(^(x')W^(x))) d (A77) = J V[A]U Av [A\ Tt(w\[A{x')]W a [A(x))) (58) 

where CZa^A] is a functional which specifies the probability to have a given field configuration in the target qq dipole. 
When A-q ~ 0, there are no extra gluon radiation and the amplitude provides the initial condition to the JIMWLK 
equation Q that describes how C/aiJ-A], and thus also the target average 1)58(1 ■ varies as one increases the rapidity 
interval Aij. A lot of work is employed trying to solve the JIMWLK equation, at least numerically but this is far 
beyond the scope of this paper. Providing the rapidity interval A77 is not too high, the JIMWLK equation reduces to 
the BFKL equation whose solution has been known for years. A first solution valid beyond BFKL has been found in 
|28j but it is only valid as long as the color fields in the two dipoles are weak. 

Anyway, as shown by formula (|5T)l . the forward-jet measurement is a suitable observable to study the gg dipole-qq 
dipole scattering and saturation effects purely due to quantum evolution. Forward-jets at HERA seem to be the good 
candidate to perform the present studies |l^ . 



C. Mueller-Navelet jets at hadron colliders 

In the future, Mueller-Navelet jets 0] at LHC could be decisive while at the Tevatron the rapidity range is probably 
not large enough 0, 0] . Mueller-Navelet jets are processes in which a proton strongly interacts with another proton 
or antiproton and where a jet with a hard transverse momemtum is detected in each of the two forward directions, 
see Fig. 3b. Such processes can also be computed from formula l|52[): the target is itself a hadron emitting a forward 
jet. The two cuts on the transverse momenta of the jets will be denoted Et x and Et 2 and Aij is now the rapidity 
interval between the two jets. One has 

9p{xj 1 ^ 2 )g P {xj 2 ,ii 2 ) I d 2 rd 2 z (j)j(r,E Tl )<f>j(z,E T2 )<j {gg){gg) (r 7 z 7 ATi) (59) 



dx j l dx 



and now the cross-section V( a g)( a g) for the scattering of two gg dipoles is the quantity driving the evolution with the 
rapidity. 1(57(1 and (|59() are generalizations of the known formulae for forward jets ^jjj and Mueller-Navelet jets 0] 
that included only linear BFKL evolution. Full quantum evolution is included in ((57(1 and ((59(1 via the gg dipole-gg 
dipole cross-section. That makes those two measurements suitable observables for studying saturation due to evolution 
and unitarization of the dipole-dipole scattering, as well as for testing phenomenological models [l3|. 



p (p) 



a) 



b) 



FIG. 3: a) Forward-jet production in a proton-lepton collison (at HERA), b) Mueller-Navelet jet production in a proton-proton 
(at LHC) or proton-antiproton (at Tevatron) collision. Et, Et 1 , and Et 2 ■ cuts on the transverse momenta of the jets; x j, xj 1 , 
and x j 2 : longitudinal momentum fraction of the jets with respect to the incident hadron; A77 : rapidity intervals between the 
hard probes. 



We have derived inclusive and diffractive forward-gluon production in the scattering of a qq dipole off a unspecified 
target in the high-energy eikonal approximation. We worked in a framework in which the incident qq pair is represented 
by its light-cone wavefunction and in which the interaction with the target is described by Wilson lines that shift 
the components of the incident wavefunction with eikonal phases. We have also assumed that the emitted gluon was 
soft so that our results are valid at leading-logarithmic accuracy for this measured gluon. The results are given by 
formulae l|18(l and l|27(l for the inclusive and diffractive cross-sections respectively. All multiple gluon exchanges with 
the target are included. The measured gluon being always the most-forward gluon with respect to the incident qq 
dipole, quantum evolution is contained in the target-averaged quantities. For the inclusive cross-section, only one 
such quantity is involved: the forward scattering amplitude of a colorless pair of gluons on the target fH^I . In the 
diffractive case however, more complicated qqg amplitudes play also a role l|2()l) . 

In the inclusive case, without specifying the target, we were able to simplify the result (|18fl and to write the cross- 
section in the factorized form (|40|l . It expresses the forward-jet cross-section as the convolution between an effective 
dipole distribution (|41|l and the gg dipole-target cross-section. The dipole factorization (|4C]|) for a forward-jet emission 
is the equivalent of the usual dipole factorization for a virtual photon: this is a dipole formalism for a hadronic probe. 

As a simple application of this result, we described forward jets at hadronic colliders. We extended the result to the 
case of an incident hadron by including linear evolution before the emission of the measured gluon jet and then by using 
the collinear factorization properties of QCD. Formula (j55l) is thus is very simple dipole-factorized cross-section that 
describes forward jets at HERA, Tevatron or LHC, depending on the target. We finally gave the expression for those 
forward-jet cross-section (|57j) and Mueller-Navelet jet cross-section l|59fl. They are written in terms of the gg dipole-qq 
dipole or gg dipole-gg dipole cross-sections which are the quantities that determine the high-energy behavior via their 
rapidity evolution. We have therefore established a link between two observables measurable at HERA, Tevatron or 
LHC and dipole-dipole cross-sections which are the basic quantities to study to understand high-energy scattering in 
QCD, saturation and unitarization due to quantum evolution. 
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